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Abstract 

We consider an open quantum system of N not directly interacting spins 
(qubits) in contact with both local and collective thermal environments. The 
qubit-environment interactions are energy conserving. We trace out the variables 
of the thermal environments and N — 2 qubits to obtain the time-dependent re- 
duced density matrix for two arbitrary qubits. We numerically simulate the re- 
duced dynamics and the creation of entanglement (concurrence) as a function of 
the parameters of the thermal environments and the number of qubits, TV. Our 
results demonstrate that the two-qubit entanglement generally decreases as N in- 
creases. We show analytically that in the limit N -» oo, no entanglement can be 
created. This indicates that collective thermal environments cannot create two- 
qubit entanglement when many qubits are located within a region of the size of the 
environment coherence length. We discuss possible applications of our approach 
to the development of a new quantum characterization of noisy environments. 

1 Introduction 

In open many-body systems, such as solid-state and biological ones, macroscopic quan- 
tum behavior reveals itself in many ways. Often the quantitative parameter used to 
measure "quantumness" (possibly of macroscopic order) is entanglement. The presence 
of entanglement implies that the wave function (or the reduced density matrix) cannot 
be represented as a product of the corresponding objects for the individual qubits. It 
is important to note that to produce and to measure entanglement in such systems, 
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one does not necessarily need to know much detail about the system, possibly not even 
its Hamiltonian [1]. The questions then are: How useful is entanglement as a measure 
of quatumness, what can it add to our knowledge of system properties and behavior, 
and how can it be utilized? Indeed, just knowing that the system is entangled (knowl- 
edge of complicated quantum behavior) is not sufficient to imply that its quantum 
properties are useful for specific applications. Fortunately, however, in some situa- 
tions entanglement provides very useful properties, including additional exponential 
resources for quantum computation [2] and a possible enhancement of photosynthesis 
in bio-systems [3]. (See also references therein.) 

Entanglement could be produced by direct interaction between the qubits. This 
interaction should be of a "conditional" nature, mixing an initial product state (dis- 
entangled) in such a way that the final state becomes correlated in a quantum way. 
Many aspects of entanglement creation are widely discussed in the literature. (See, for 
example, [4-9] and references therein.) 

More recently, interest appeared in the possibility to create entanglement in the ab- 
sence of direct interactions between qubits (or when the latter are very small). Entan- 
glement can then still be created merely by the indirect interaction of non-interacting 
qubits through a collective thermal bath. In [4] this situation was considered for a 
model of two non-interacting spins 1 /2 (qubits) interacting only with a collective ther- 
mal bosonic environment. It was demonstrated numerically in [4] that for some initially 
unentangled two-qubit states, and under some conditions on the the thermal bath, 
measurable entanglement between the two qubits is created for intermediate times. 
The model of [4] is energy-conserving, ignoring relaxation processes for the qubits, 
and including only the effects of decoherence. In [9] these results were extended to 
a more general model having (i) both local and collective thermal environments (at 
the same temperature) and (ii) energy conserving and energy exchange interactions 
between qubits and their environments. The conditions for entanglement creation were 
discussed and analyzed numerically in [9]. It was concluded that, in spite of the com- 
petition between the local thermal environments (which destroy entanglement) and the 
collective thermal environment (tending to create entanglement), the creation of mea- 
surable entanglement can be realized for some finite times. In both papers, [4] and [9] , 
only two qubits are analyzed, hence no direct connection to many-body systems was 
made. As was recently shown in [10], the presence of a large number of indirectly inter- 
acting qubits interacting only with their common collective thermal environment could 
significantly modify the effective single-qubit characteristics including their relaxation 
and decoherence rates. 

In the present paper we propose a new approach to the analysis of thermal envi- 
ronments, namely, the analysis of their ability to create quantum entanglement. For 
example, in biological systems, noise produced by a thermal environment can usually 
be described using standard noise characteristics, including space and time correlation 
functions (with the corresponding correlation lengths), the spectral density of noise, 
and high order correlations (if needed). These characteristics allow one to introduce, 
estimate, and measure important parameters of the biological system, such as relax- 
ation and decoherence (dephasing) times, kinetic and transport coefficients and others. 
We propose to take the next step in this direction: to develop a new method for char- 
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acterizing noisy environments created by biological species at room temperature - a 
quantum characterization of noisy environment. 

The basic parameter measured in our method is the two-qubit (two-spin) concur- 
rence describing the level of two-qubit entanglement created by noisy environments 
at room temperature. Experimentally, the method can be implemented using the 
well-developed technique of liquid state nuclear magnetic resonance (NMR) quantum 
computers, using an ensemble of two-qubit systems. The qubits are represented by 
nuclear spins in molecules which are diluted in a liquid. The method we propose will 
allow one to determine the level of entanglement of two spins by a particular noisy 
environment. Different types of noisy environments can be distinguished using this 
quantum analysis. This approach will allow one to develop new quantum technology 
to characterize the noisy environments created in biological systems, by measuring the 
level of entanglement (concurrence) the environments generate between two nuclear 
spins. 

Suppose that an ensemble of two-spin molecules is diluted in a solution filled with 
microorganisms (e.g. bacteria or viruses). We single out two nuclear spins, which can 
be studied by conventional NMR techniques at room temperature. So, in our ensemble 
each two-spin molecule is attached to a single bacterium (virus). We assume that 
these chosen molecules will attach to bacteria so that the two spins of the attached 
molecule interact through the noisy environments (local and collective) produced by 
the bacteria. One considers three such samples of diluted two-spin molecules attached 
to bacteria. Three samples are necessary for implementing a spatial labeling [2]. (We 
choose spatial labeling for definiteness, but temporal labeling with a single sample is 
possible as well.) External electromagnetic pulses can selectively excite each sample. 
The signals from the three samples are added to obtain their average. In the method 
of spatial labeling, the spin systems in all samples are initially in equilibrium at room 
temperature. Then a preliminary sequence of radio-frequency (r/) pulses is applied to 
two of the three samples in order to change their mixed state by permuting diagonal 
density matrix elements. Then the evolution of the three-sample average NMR signal 
is exactly the same as if the two spins were initially in their ground states. One can 
reasonably doubt the possibility of entangling spins which begin their evolution from a 
state of thermal equilibrium with the noisy environment. However, labeling (averaging) 
reveals the evolution of the NMR signal, which corresponds to the initial ground states 
of the spins. One can consider this spatial labeling as the simulation of the NMR 
signals, that would be obtained if the spins were initially in their ground state. Thus, 
in this case one may consider an ensemble of two-spin systems, which start from the 
ground state. 

After the initial preparation of the samples (permutation of the density matrix 
elements), all three samples are subjected to the same sequences of the rf pulses. First, 
one applies initial pulses which change the mutual orientation of the spins. Then, 
the noisy environment produced by bacteria can generate entanglement. The level 
of entanglement, which may be negligible for one kind of bacteria and significant for 
another kind, is determined by measuring the concurrence. To find the latter, one 
applies, to all three samples, electromagnetic pulses to carry out a quantum state 
tomography protocol [2]. The state tomography protocol allows one to recover the 
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spin density matrix. (More precisely, in order to implement the state tomography 
protocol, one repeats the whole procedure nine times with different sequences of state 
tomography pulses.) Once the spin density matrix is recovered, one is able to compute 
the concurrence by measuring the entanglement created by the noisy environment. This 
determines whether the noise produced by bacteria creates quantum entanglement in 
a two-spin system. We assume that the scalar J-J interaction between the two nuclear 
spins mediated by the chemical bonds in a molecule does not produce entanglement. To 
achieve the latter case, the mutual orientation of the spins created by the initial pulses 
should be chosen in such a way that the J-J interaction does not produce entanglement. 

In order to be feasible this approach requires the following: (1) The interaction with 
the noisy environment must produce significant entanglement during a time shorter 
than the relaxation and decoherence times. (2) The correlation length of the noisy 
environment must be greater than the distance between the two spins in the molecule. 
(3) The relaxation and decoherence times in the two-spin system must exceed the time 
required to apply the pulse sequence in the quantum tomography measurement. 

Condition (1) imposes a requirement on the effective constant of the spin-environment 
interaction. Condition (2) imposes a requirement on the correlation radius of the spin- 
environment interaction. Condition (3) is satisfied in liquid NMR [2]. 

As a concrete example one can consider molecules of chloroform [2] . In chlorophorm, 
the frequencies of the hydrogen and carbon nuclear spins are well separated and can 
be addressed independently. The dipole-dipole interaction between the spins in liquids 
is suppressed due to the motion of the molecules. The J-J scalar coupling is mediated 
by chemical bonds with a corresponding frequency of about 200 Hz. In this example, 
an initial non-equilibrium state should be chosen so that the J-J coupling does not 
create its own entanglement with a significant value of concurrence. The relaxation 
and dephasing times for the proton spin are 18s and 7s, and for the carbon spin 25s 
and 0.3s. These times are long enough to implement the two-spin tomography. 

We note that any quantum information processor could be used for the study of 
entanglement induced by a noisy environment. One could consider, for example, two 
ions in an ion trap quantum computer [2] or two superconducting qubits [11]. Our ex- 
ample with the NMR quantum information processor is associated with the possibility 
to study entanglement induced by biological objects and has as general goals the cre- 
ation of the novel devices exploiting quantum effects at room temperature. The NMR 
quantum information processor is, probably, the only recently available technique of 
this kind. 

In the present paper, we consider a model of ./V not directly interacting spins 1/2 
(qubits) placed in a constant effective magnetic field (oriented in the z-direction). The 
qubits interact with both local and collective thermal environments (all at the same 
temperature). The collective interaction introduces an indirect qubit interaction. In 
the total density matrix of all qubits and environments, we trace over the variables of 
the environments and N — 2 qubits. This gives us the time-dependent reduced density 
matrix for two arbitrary qubits. In the ^-representation it is represented by a time- 
dependent 4x4 matrix. It is important to notice that the matrix elements, [pt\ n ,mi 
n,m = 1,...,4, depend not only on the parameters of the thermal environments but 
also on the total number of qubits, N. We study numerically the concurrence C(t) 
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of the reduced two-qubit density matrix, and its dependence on the parameters of the 
system and on N. To realize and study this situation in an experiment, one must 
have access to the two selected qubits (such as their particular frequencies) , in order to 
manipulate them and prepare the initial state. Our main result is that the amplitude of 
concurrence, C max (t), generally decreases as N increases. This means that one should 
not expect that the collective thermal environment can create by itself measurable 
entanglement even of two-qubits, in the presence of many other qubits within a range 
of the collective environment coherence length. 

1.1 Outline of main results 

The initial state of the entire system is disentangled, a product state in which each of 
the N spins is in a state pj, j = 1, . . . , N, all local reservoir states are thermal and so 
is that of the collective reservoir (at a fixed common temperature x ). 

Analytic results. 

• Explicit dynamics. As the spins interact with the reservoirs via energy-conserving 
couplings only, the reduced two-spin dynamics can be calculated explicitly, see Propo- 
sition 2.1 below. Consequences of the energy conservation are that populations, i.e., 
the diagonal density matrix elements, are time- independent, and that the off-diagonal 
elements evolve independently. As an example, we discuss here the dynamics of the 
(1,2) matrix element, 

[Pt]i2 = [po]i2 j* s ® e -*!m)-»lrc{t) p N{ty (L1) 

The other matrix elements have similar behavior. Each factor on the r.h.s. has an 
interpretation: 

- [po]i2 is the initial condition of the matrix element in question. None of the other 
initial matrix elements are involved (energy conserving coupling); 

- e 1 ^ 2 * is the uncoupled dynamics (no interaction with environments); 

- e 1KcS ^ is a dephasing factor with a time-dependent phase S(t) < becoming 
linear for large t (for the considered infra-red behavior \k\ 1 / 2 of the coupling constants 
in three dimensions, see (2.9)); it represents a "Lamb shift" contribution to the real 
part of the effective energy; this term is generated by the collective reservoir, but it is 
independent of the presence of the N — 2 traced-out spins (the term would be the same 
if only two spins were coupled to the reservoirs); 

- e~ K e r (( t )- K c r c{t) j g a decaying factor with time-dependent decay rates, T(t) > 0, 
becoming linear for large t (see (2.10)) both the local and collective reservoirs con- 
tribute; however, the term is independent of the N — 2 traced-out spins (again, it 
would be the same if only two spins were coupled to the reservoirs); 

- P/v(i) is a product of N — 2 oscillating terms encoding the effect of all the traced- 
out spins (see (2.8)). It is important to notice that -Pjv(i) only depends on the diagonal 
density matrix elements of the initial states of the N—2 traced-out spins. 2 Consequently, 

1 A generalization to a non-equilibrium situation where each local and the collective reservoir have 
different individual temperatures is immediate. 

2 This is so since the dynamics is energy-conserving, and tracing out any of the spins involves only 
the diagonal of the initial (time zero) density matrix. (See also Remark 1 after Proposition 2.1.) 



5 



the two-qubit state does not depend on the initial off-diagonal density matrix elements 
of the N — 2 traced-out "background" spins. Typically, we expect those spins to be 
initially in (close to) equilibrium, corresponding to vanishing off-diagonals. 

Some general properties of Pjv(i) can be explained easily for the case in which all 
N — 2 spins are initially in the high temperature equilibrium state ( | +) + | — ) ) . Then 

P/v(i) = [cosf^c S(t))] N ~ 2 an d its magnitude oscillates between zero and one. Due to 
the large power (N — 2), the peaks of the function |Pjv(i)|, centered at the discrete 
times t p satisfying S(t p ) G -^-Z, are of very narrow width 0(1 /(x^\/iV)) for large N. 
Consequently, in the limit N — > oo, with x c held fixed, |P/v(t)| is zero for all t except for 
t = t p , where |Pjv(i p )| = 1. But the density matrix becomes very simple if P/v(t) = 0, 
because many entries vanish (c.f. Proposition (2.1)) and the corresponding concurrence 
is zero. It follows that in the large N limit, concurrence is zero for all times (except 
possibly for some isolated instances, t p ). 

• N -dependent scaling of the interaction x c . The above analysis suggests that 
one cannot generate two-spin entanglement for N large at fixed interaction strength 
x c . However, the width of the peaked function PN(t) which is of order 0(l/(x^y/N)) 
becomes appreciable if x c > N 1 / 4 . Hence we consider a iV-dependent scaling of the 
coupling, replacing x c by K c jN n ', for some r/ > 0. According to the above discussion, 
the borderline case is r/ = 1/4. 

Starting from the explicit expressions (Proposition 2.1) and using the scaling >c c /N v , 
we calculate the limit N — > oo of p t , for t € R fixed. The analytic expressions we 
obtain for < r/ < 1/4 and w > 1/4 show that the limiting dynamics does not create 
entanglement, for any time t. While we are able to obtain explicit expressions for 
concurrence in the regime of N — > oo, we are not so for iV finite. 3 However, since no 
entanglement is generated in the limiting case, iV — > oo, but we know entanglement is 
created for iV = 2 (see e.g. [4,9]), we expect that entanglement creation decays with 
increasing N. We study this decay numerically. 

Numerical results. 

We introduce f c , the highest frequency at which spin-reservoir interactions occur 
and call it the cutoff frequency. In the simulations, we take v c of the order of the 
thermal frequency vt = ksT/h. In the infra-red regime, our coupling is proportional 
to ^/\k\ (see after 2.10). 

• For N = 2, concurrence creation is maximal if both spins start out in their 
high-temperature state ^(|+) + | _ )), see Fig.l. Consequently, in the subsequent 
simulations, we take initial states of the two not traced-out qubits very close to this 
state, and we take the diagonals of the initial states of the N — 2 traced-out quibts to 
be constant 1/2 (remember that the off-diagonals of these qubits do not influence the 
dynamics at all). In Fig.4 we modify the initial state of the two not traced-out qubits 
and check that maximal concurrence is indeed obtained when both qubits are in the 
above state, even for large N. 

• For general N, entanglement evolves according to a rescaled time t h-> x^f c t, see 
Fig. 2. This figure shows that a reduction of x c diminishes the created concurrence in 

3 The reduced density matrix is given explicitly for all N and all t, but calculating from it explicitly 
the concurrence is more difficult. 
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a moderate way. For instance decreasing x c by a factor 10 only decreases concurrence 
by less than 1/3. 

• In Fig. 3 we show that the maximum of created concurrence decays with increasing 
N. For intermediate values of N (with the current parameters N ~ 10 — 150) the 
decrease is exponential, for smaller and larger values of N, it is superexponential. 

• In the same Fig.3c, we study the dependence of the maximal time, r c , (before 
recurrence) at which the concurrence is not zero. We have found that this time decays 
exponentially in the number of spins, N, for sufficiently large N. 

• Results on the rescaled model x c i— > x c /N v are shown in Fig. 5. We find a decrease 
of maximal concurrence with increasing N for all rj > 0. The critical value, rj = 1/4 (see 
analytic results above) divides the concurrence decay into two regimes. In the range, 
rj > 1/4, the maximal concurrence decreases exponentially in N, for intermediate values 
of -/V (between 10 and 180), with a universal decay rate (i.e., not depending on rj). For 
77 < 1/4 the decay is superexponential and varies with r\. We conclude that no scaling 
x c i y Xc/N^ can compensate the decay of created concurrence for large N. 



2 Model and reduced density matrix 

The full Hamiltonian of the N noninteracting spins 1/2 coupled by energy conserving 
interactions to local and collective bosonic heat reservoirs is given by 

N N 

H = -hJ2"nS* + J2 H *n + HR (2.1) 

n=l n=l 

N N 

+ *nS Z n ® (f> c (/ c ) + U nSn ® 4>n{fn)- (2-2) 
n=l n=l 

Below we use dimensionless variables and parameters. To do so, we introduce a 
characteristic frequency, uq, typically of the order of spin transition frequency. The 
total Hamiltonian, energies of spin states, and temperature are measured in units fvjjQ. 
The frequencies of spins, oj n > 0, bosonic excitations, uj(k) = c\k\ (where c is the 
speed of light), the wave vectors of bosonic excitations are normalized by ojo/c, and all 
constants of interactions are measured in units of ujq . The dimensionless time is defined 
as uiQt. 

In (2.1), (2.2), uj n > is the frequency of spin n, 



(2.3) 



b ~ 2 [ -1 

and S* denotes the S z of spin n. Hr is the Hamiltonian of the bosonic collective 
reservoir, 

H R = [ \k\a*(k)a(k)d 3 k, (2.4) 

and H Rn is that same Hamiltonian for the n-th individual reservoir. For a square- 
integrable form factor h(k), k G M 3 , 4>{h) is given by 

<t>{h) = ^=l {h(k)a*(k) + h(k)*a(k)}d 3 k. (2.5) 
V2 7r3 
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The real numbers, x n and v n , are coupling constants, measuring the strengths of the 
energy-conserving collective coupling and the energy-conserving local coupling, respec- 
tively. 

Since the spins interact with the reservoirs only through energy-conserving channels, 
this model is exactly solvable. For simplicity of exposition, we take 

x n = k c for all n (collective) 

and 

v n = v i for all n (local). 

We also take, for simplicity, all local form factors equal (ft) and all collective ones also 
(/,)• 

Fix any pair of spins, and (re-)label their frequencies by uj\ and u)2, see (2.1). We 
write the reduced density matrix, pt, of the two fixed spins as a 4 x 4 matrix [pt]ij in 
the ordered energy basis 

$i = (fi <8> <pi, $2 = <pi <8> <p-i, $3 = <P-i <8> <fi, $4 = <P-i <8> (f-i, 4 (2.6) 

where S z ip± 1 = ±§<p±i- For instance, [pt\2,A = ($2,Pt$4)- 

The initial state of the spins is a product state of the form ps lt o <8> • • • <8> Ps N ,o, where 



PS j,0 



Pj v j 
L V* 1 - Pj 



(2.7) 



with < pj < 1 and \vj\ 2 < pj(l—pj). The upper bound on the off-diagonal guarantees 
that the eigenvalues of pSj,o are non- negative. 
We introduce the quantities 

N 

P N (t) = J] [ Pj + (1 - Pj ) e-^^W] , (2.8) 



J=3 
1 

2 



sm = 4 / \m\ 2 lklt -Z m) i% (2.9) 



r,, c (t) = / |/ t , e (m 2 coth(g|t|/2) ^ ^l t/2> d 3 fc (2.10) 

The integrals in (2.9) and (2.10) are made to converge introducing a suitable cut- 
off wavenumber, \k c \, or cut-off frequency, v c = \k c \/2ir. (Here we use dimensionless 
units.) For instance, for numerical simulations, we choose as form factor the function 
fc(k) = vW X|Jfc|<|fech where X|fe|<|fe c | = 1 if \ k \ ^ \ k c\ and X|fe|<|fe c | = otherwise. 

We also define -Pjv(t) to be the same as Pjv(*)> bvrt with x% replaced by 2x^. With 
this notation, we have the following result. 



4 Another equivalent notation is: $1 = | + + >, $ 2 = | H — >, $3 = | - + >, $4 = >• 
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Proposition 2.1 (Explicit dynamics of the reduced density matrix.) The evo- 
lution of the density matrix is given by 



[Pt] 12 = 


[Po]l2 


e i^s(t) e -x|r £ (t)-x 2 r c (t) p N ( t ^ 


(2.11) 


[ft] 13 = 


[ft] 13 


jwt e ix?s (t) e -x?r c (t) Pjv(t)) 


(2.12) 


[pt]u = 


[ft] 14 


e iK+^)< e -2x|r,(t)-4x 2 r c (t) p^f^ 


(2.13) 


[ft] 23 = 


[ft] 23 


e i(u;i-u;2)« e -2x|r £ (<) ) 


(2.14) 


[ft] 24 = 


[ft] 24 


junt e -i^s(t) e -^r e (t)-^r c (t) p N u\ 


(2.15) 


[ft] 34 = 


[ft)] 34 


e iw 2 t e -ix 2 c s(t) e -^r e (t)-x 2 c r c (t) p N ( t ^ 


(2.16) 



and the populations are constant, [pt]jj = [po]jj> for j = 1, . . . , 4 and ( G R. 

The proof of this proposition is a rather simple calculation. One can proceed as 
in [12] (proof of Proposition 7.4). 

Remarks. 1. The effect of spins 3, ...,7V is contained entirely in the factors 
P/v(t) and -P/v(t). They only depend on the initial populations pj, j = 3,...,N (see 
(2.8)), but not on the off-diagonals, Vj. This is explained by the fact that when tracing 
over a single spin, j > 3, we perform the operation Tr spin j Ups j: oV, where U, V are 
operators commuting with Sf (energy conserving interactions only!). Clearly the latter 
trace only involves the diagonal of ps $■ 

2. The oscillatory phases, e lujt , in (2.11)-(2.16) represent the free, uncoupled dy- 
namics of the spins. Consider the modified two-spin density matrix 

p> — gitf-wiSf-G^Sf) pt g-itC-wiSf-G^Sf^ (2.17) 

("interaction picture" dynamics of pt). Because p' t and pt are related by conjugation of 
a unitary operator of the product form, e 1 *( _a,1,s 'i) <g) e 1 * (_W2 ' s 2) ) the concurrences of p t 
and p' t are the same. In other words, when examining concurrence of pt, we may use 
formulas (2.11)-(2.16) with uj\ = u>2 = 0. 



2.1 Concurrence 



Recall that the concurrence of the reduced density matrix, pt, is unchanged when we 
pass to the interaction picture p' t (see the remark explaining (2.17)). In the basis (2.6), 
the evolution of p' t , (2.17), is given by (2.11)-(2.16) with oj\ = U2 = 0, and where the 
initial condition is (p' = po) 



[Po 



P\P2 PlV2 
P\V* 2 ft(l-P2) 
V{p2 V\V2 

P2) 



V\P2 V\V2 
V\V\ V\(\~P2) 
(1-Pl)p 2 (1-Pl)f2 

(i-piK (i-ft)(i-p 2 ) 



2.1.1 Variation of iV-dependence 

For homogeneous initial conditions, pj = p for j = 3, . . . , TV, we have 

P N (t) = [ pe ^ s ^ + (1 - p)e~ i ^ s(t Y~ 2 - 



(2.18) 



(2.19) 
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Unless p = 0,1, |Pjv(i)| oscillates in t between its minimum value \2p — 1\ N 2 (when 
cos(2x 2 S(t)) = -1) and its maximum value 1 (when cos(2x 2 S(t)) = 1). 5 The width of 
the oscillations becomes very narrow with increasing N. In the limit of large N, Piv(t) 
is zero for all times, except for the discrete set of t G R satisfying cos(2x 2 S(t)) = 1, in 
which case |P/v(t)| = 1. 

This implies that for large N, all off-diagonal density matrix elements of pi vanish 
with the exception of [pt\2?> = e 1<wl ~ W2 )*e~ 2x ? r ^*) (and \pt\32, of course) for almost all 
values of t. This suppression of off-diagonals comes from the large number of particles 
and is mediated through the collective energy-conserving interaction. (For x c = we 
have P N (t) = 1.) 

In order to try to have a non-trivial dynamics for large N, one may scale the 
collective conserving coupling constant as 

k c — > — — , some 77 > 0. 

Then (2.19) becomes 

P N (t) = \p^mi^> + (1 _ p)e -^ c s(t)/N^ ]N ~ 2 . (2.20) 
An expansion in large N yields 

P N (t) = e -i* 2 cS(t)ll~2p]Ni-^ e -2^(t)N^[p(l- P )+0{N-^)] 

oo, 




e -i^S(t)[i-2p]ivV2 e -2*isHt)p(i-p)N^ o < t? < 1/4 (and p ± 0, 1) 

e -i* 2 S(t)[l-2p]iV 1 - 2 ') 1/4 < 7/ 

(2-21) 

Remarks. 1. By replacing, in these limits, x 2 by 2x 2 , we obtain the corresponding 
limits for Pjv(i). 

2. For < 77 < 1/4, Piv(t) vanishes as A — >■ 00. 

3. The rapid oscillating phases disappear if p = 1/2 (any 7/ > 0) or 7/ = 1/2 (any 

p). 

2.1.2 Asymptotic concurrence (A — >■ 00) 

• < 77 < 1/4 : The reduced two-spin density matrix (in the interaction picture) at 
time t is 

P1P2 

Pi(l-P2) vi^e -2 *? 1 ^ 
vlv 2 e- 2x e r tW (l-pi)l>2 
(I-P1XI-P2) 



[ft] 



(2.22) 



from which we obtain the concurrence 

C(p t ) = max {0, -2[Vpi(1 1 Pi)p 2 (l 1 P2) - |«i| Me" 2 ^^]} = 0. (2.23) 
5 We have |pe i,f ' s(t) + (1 - p) e - i ^ s( ' ) | 2 = p 2 + 2p(l - p) cos(2x 2 S(t)) + (1 - p) 2 . 
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Remember that |t>j| 2 < pj(l — pj)- This shows that the N — > oo asymptotic dynamics 
cannnot create entanglement at any time. 

• r] > 1/4 : Call the r.h.s. of (2.21) Poo(t) (a quantity still depending on TV unless 
rj = 1/2). By replacing x 2 by 2x 2 in (2.21) we obtain the limit of Pjv(i), which we 
call Poo(t). F° r V > 1/4 we have the relation Poc(t) = [Poo(t)] 2 . The reduced two-spin 
density matrix (in the interaction picture) at time t is, for N — > oo 



[Pt] 



PlP2 

PivZDtPnity 
vfaDtPooity 

vMDjiP^t) 2 }* 



p^DePooit) v^D^Pooit) 
p 1 (l-p 2 ) viv^Dj 
v\v 2 Df (1 - p ± )p 2 

vt(i- P2 )D e p 0O (ty (l-piMDtPooity 



(2.24) 

v^DjP^tf 
Vl (l - p 2 )D e P 00 (t) 

(1-P1)V2D/Poo(t) 
(1-Pl)(l-P2) 



where 



D e = D e {t) = e-^ l \ Poo{t) = e 
The density matrix (2.24) is of the product form 



-ix2S(t)[l-2p]AT 1 - 2 " 



Pl «l^(t)Poo(t) 

^^(OPooW* 1-pi 



P2 v 2 D e (t)P OD (t) 

V*D e (t) P^t)* 1-P2 



This shows that the iV — >■ oo asymptotic dynamics is factorizable and cannnot create 
entanglement at any time. 



3 Numerical Results 

• Let us first consider the case of two spins only, N = 2. In (2.11)-(2.16) we put 
for simplicity = T c = T and regard T and S as two independent parameters. Taking 
both spins initially in the same state given by p, v, see (2.7), we examine the maximal 
concurrence, as a function of S and T, for arbitrary fixed values of p and v. 

We find that for fixed p, v, the maximal concurrence is given at S = 7r/2, T = 0. 
Having such values fixed and plotting the concurrence as a function of p, v, the maximal 
concurrence is realized when p = v = 1/2, see Fig. 1, where a plot of the concurrence 
as a function of p, v is shown. Maximal generation of concurrence is thus obtained 
starting from pure state initial conditions -75G+) + I - )) f° r each spin. 

• Let us now consider the case of N spins. For concreteness we choose p = 1/2 
for all spins (the traced-out ones and the two not traced-out ones). For the two not 
traced-out spins we take off-diagonals v = 0.48. (Then p is close to v which favors larger 
entanglement creation). Recall that the dynamics is independent of the off-diagonals of 
the N — 2 traced-out spins (i.e., we do not have to specify the v of the N — 2 traced-out 
spins). 

As mentioned after (2.10), we choose the form factor, f c (k) = \f\k\ X\k\<i-Kv c i with 
the cut-off frequency equal to the thermal frequency, v c = vt = ksT/h, at room 
temperature, T = 300 K. 

In Fig. 2 we investigate the effects of an increase in the coupling parameter, x c . 
The first effect is a time-shift for the concurrence evolution, described by a scaling, 



11 



t —7- x^v c t, see Fig. 2a, where v c is the cut-off frequency. The second effect is a 
reduction of the maximal concurrence in a smooth way, see Fig. 2b. As one can see, 
the effective decrease in amplitude for not too strong coupling strength, x c , is relatively 
small. For instance, changing x c for TV = 2 by one order of magnitude from 0.04 to 
0.4, changes the amplitude by only 27%. The percentage change is almost the same 
for larger TV values, see Fig. 2b. 

In Fig. 3 we show that the creation of concurrence decreases with the number of 
spins. In graph a) we plot the concurrence as a function of (rescaled) time for various 
values of TV = 2, ... , 32. As one can see, the same time rescaling is also valid for TV > 2. 
Moreover, the maximum concurrence created, C max , reported in graph b), decreases 
exponentially in TV in the range 10 < TV < 150 and faster than exponentially outside 
this range. For larger TV, the concurrence decays superexponentially in TV. For TV 
exceeding 200, the concurrence becomes too small to be significant (of the order 10~ 4 ). 

It is also interesting to note that the graph of concurrence shows collapses and 
revivals and that the revival times for TV > 2 are always less than the revival time 
for TV = 2. It is also interesting to consider how the collapse time, r c , defined as 
the first time at which concurrence drops abruptly to zero, depends on the number of 
spins, TV". This study has been reported in Fig. 3c and shows that the rescaled collapse 
time decays exponentially with the number of spins: r c = M^,v c t c ~ exp(— aN), where 
a = 0.0838 ± 0.0002. (See dashed line in Fig. 3c.) 

• One can also vary the initial conditions for the spins by choosing independent 
pi t 2 and V12, while all other spins have the same value pj = 1/2, j = 3, . . . , TV (their 
off-diagonals Vj do not influence the dynamics at all). 

In order to simplify the problem, we also set p\ = v\ and p2 = V2 and consider the 
maximal concurrence as a function of two independent parameters, p\ and p2, only. 

An example of the 3D plot obtained is reported in Fig. 4 : the maximal concurrence 
is realized for p\ = p2 = v\ = v% = 1/2, independently of TV (in the picture TV = 40, 
but similar plots are obtained for other values of TV) . 

• The numerical analysis of the rescaled model with x c replaced by x c /N v shows 
that the concurrence is always a decreasing function of TV and that the maximum of 
the created concurrence is a universal function of the number of spins TV, independent 
of 7] for T] > 1/4. 

Results are shown in Fig. 5, where the dashed line is the best exponential fit 
exp(— aN), with a = 0.0177 ± 0.0003, the best fitting value, for the cases r\ > 1/4. 
The same figure shows that when rj < 1/4 the decay is superexponential and no uni- 
versality occurs. 

This suggests that no power law scaling with TV of the coupling strength can com- 
pensate the rapid decay of concurrence with the number of spins. 
Qualitatively similar results, not reported here, can be obtained by changing the ratio 
between the thermal and cut-frequency in the range (0.5,4). 

4 Conclusion 

We present a new way to characterize a noise source by analyzing its ability to create 
entanglement between two arbitrary qubits in the TV-qubit open system. We have 
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discussed an application of this approach to the analysis of noise produced by bacteria 
or viruses. As a first step we consider the dependence of the concurrence on the 
number of qubits. We show that concurrence quickly decays with increasing number 
of surrounding qubits. It follows that for implementing our approach, one has to 
use a small number of qubits collectively interacting with the thermal environment, 
preferably only two qubits. 
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N= 2 spin r = 0, S = tt/2 




o 



Figure 1: Maximal concurrence as a function of p and v, for fixed T(t) = ans S(t) — ir/2. 
Here, N — 2 spins is considered. 
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Figure 2: a) Concurrence as a function of the rescaled time, for different >c c and fixed e = 
VcI^t = 1) being respectively the cut-off and thermal frequency (yi — ksT/h with T = 300-ftT) 
and the form factor f c (k) = Other data are p = 0.5, v = 0.48, 2V = 4. b) Plot of 

the concurrence at the peak (obtained from a)) as a function of the coupling strength, k c , for 
different N values as indicated in the legend. Other values are the same as in a) . 




Figure 3: a) Concurrence as a function of the rescaled time, for fixed e = v c /vt — 1, being 
respectively the cut-off and thermal frequency [vt — ksT/h with T — 300^) and the form 
factor, f c {k) = \/\k\- Other data are p — 0.5, v = 0.48, xr c = 0.05. b) Plot of the concurrence 
at the peak (obtained from a)) as a function of the number of spins N. c) Plot of the collapse 
time, r c , as a function of the number of spins, N. The dashed line is the best fit exponential 
exp(-aA0 w ith a = 0.0838 ± 0.0002. 
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Figure 4: a) Maximal concurrence as a function of the independent parameters, p\ = v\ 
and P2 = V2- As one can see the maximal concurrence is realized at the external corner, i.e 
pi = p 2 = 1/2. Here is e = 1, N = 40 and p = 1/2 for all other spins. 
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Figure 5: Maximal concurrence as a function of the number of spins, for different 
power law scaling, as indicated in the legend. Here is f c (k) = \/]k~\, p = 0.5, v = 0.48, 
e = v c jvr = 1, T = 300 K, x c = 0.2. The dashed line indicates a fitting exponential 
for the cases rj > 1/4. The solid curve indicates the case, r\ = 1/4. 
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